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Abstract 

Using X-ray data, we investigate interrelations between gas density po, virial density p vlT , core radius 
R, temperature T, entropy S gas , and metal abundance Z in the core region of clusters of galaxies. First, we 
confirm that fundamental relations among po, R, and T found by Fujita and Takahara are reproduced by 
another data catalogue. Second, we find that, when clusters have two components in their surface brightness 
distribution, the inner components also satisfy the same fundamental relations on the assumption that the 
average temperature of the inner component is the same as that of the outer component. These results 
strengthen our interpretation that clusters form a two parameter family in terms of mass and p v ir! larger 
Pvir corresponds to earlier formation epoch. We argue that the inner components represent distinct dark 
matter components which collapsed ahead of the outer components. Third, we also find a tight relation 
between ,!5g as and p vlI both for the outer and inner components; S sas is smaller for larger p vlI but is larger 
than that produced through gravitational collapse alone for larger p vir . Although radiative cooling affects 
the thermal evolution, the tight relationship discovered suggests the existence of stable heating sources or 
stable energy transmission mechanisms. Finally, we find that the iron abundance at the centers of clusters 
is correlated with p v ; r for the inner components. This implies that iron produced by Type la supernovae 
has more accumulated for clusters formed earlier. We briefly discuss the implications of these findings. 
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1. Introduction 

Since clusters of galaxies are the largest virialized 
structure in the universe, they should carry the infor- 
mation on the cosmological structure formation such as 
the cosmological parameters and the initial density fluc- 
tuations. In particular, theoretical studies have shown 
that the structure of their core region should provide use- 
ful clues to these issues. Numerical simulations done by 
Navarro et al. (1997) demonstrate that the characteristic 
density of a cluster correlates to that at the formation 
redshift. By means of a modification of the extended 
Press & Schechter formalism, Salvador-Sole et al. (1998) 
show that the core radius of halos is essentially propor- 
tional to the virial radius of the clusters at the time of 
formation. So the cluster cores are expected to retain 
the properties at the formation epoch rather than to be 
completely destroyed by subsequent dynamical evolution 
after the formation. The observations of the central gas 
density po and core radius R are to be compared to these 
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predictions to give a formation redshift. 

Most extensive data sets now available to such studies 
are X-rays from hot intracluster gas. The X-ray lumi- 
nosity of clusters Lx is mainly determined by the struc- 
ture around the core region. The relation between Lx 
and X-ray temperature T has been used as a diagnos- 
tic of cosmological structure formation. However, most 
of the previous studies have not explicitly cared about 
the structure of central region, rather they avoided to 
use physical quantities there. For example, they do not 
take account of a variation of core radius for a given 
X-ray temperature. One of the reasons for this is that 
it was difficult to determine core radii because of their 
correlation with the slope of surface brightness profile 
[3, the central excess emission (so-called cooling flow), 
and the poor spatial resolutions of telescopes. However, 
ROSAT satellite makes it possible to obtain core radii 
with enough accuracy. Using the data compiled by Mohr 
et al. (1999), we have investigated the structures in 
core regions of clusters (Fujita, Takahara 1999a, hereafter 
Paper I). We find that there are fundamental relations 
among po, R, and T; clusters are distributed on a band 
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in the (log pojlog i?, logT) space forming a fundamental 
plane. Moreover, we find that the observed Lx~T rela- 
tion is the cross section of the band seen from the main 
axis of the band, and that the information along the main 
axis is almost degenerated on the Lx~T relation. By 
comparing the results with the theoretical predictions, 
we show that the main axis corresponds to a variation 
of the central virial density p v ir and thus to the cluster 
formation redshift, and that the scatter of the data along 
it implies a wide range of cluster age (Fujita, Takahara 
1999b, hereafter Paper II). In §2 of this paper, we investi- 
gate whether these results hold for Einstein data by Jones 
and Forman (1984), which is another data set giving p , 
R, and T in spite of being less accurate. 

The physical nature of the central excess components is 
another interesting issue. Recently, detailed observations 
of the central regions also become available thanks to 
the results obtained with the ASCA satellite, which has 
an imaging spectroscopic capability for 2-10 kcV range. 
Ikebe et al. (1996) point out that the central excess of 
X-ray emission seen in the Fornax cluster shows the ex- 
istence of another dark matter component which is dis- 
tinct from that representing the gravitational potential 
of the whole cluster. That is, the dark matter in this 
cluster shows a double distribution. For the Centaurus 
cluster, Fukazawa et al. (1994) and Ikebe et al. (1999) 
show that the gas attributed to the central excess emis- 
sion has at least two temperature components. Moreover, 
Xu et al. (1998) find that the central excess emission 
in A1795 is dominated by the hot component and not 
by the cold component. Although the two-temperature 
structure and/or central excess emission in many clusters 
have been interpreted in terms of cooling flows (Fabian 
1994), it is valuable to study them statistically in terms 
of a double distribution of dark matter and compare the 
results for the central emission with those for the whole 
cluster. For example, if there are indeed two compo- 
nents of dark matter distribution in clusters, hierarchi- 
cal clustering scenario predicts that the inner component 
collapsed before the outer component did according to 
the law of perturbation growing. The fundamental plane 
analysis of the inner components is done in §3. 

The discrepancy of the cluster Lx—T relation has been 
studied by many authors. While a simple scaling theory 
predicts L x oc T 2 (Kaiser 1991; Eke et al. 1996), obser- 
vations show that Lx oc T 3 (e.g. Edge, Stewart 1991; 
Markevitch 1998). One interpretation of the discrepancy 
is that the entropy of the gas in all clusters is raised at 
early times to levels comparable to those reached through 
the gravitational collapse (Kaiser 1991; Evrard, Henry 
1991). As a result, the central density has an upper 
limit; in particular, the clusters with low temperature 
are severely affected by the limit. Indeed, Ponman et 
al. (1999) find that the entropy of the intracluster gas 
is higher than can be explained by gravitational collapse 
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alone, especially for poor clusters (see also David et al. 

1996) . They argue that the excess entropy is a relic of 
the energetic winds generated by supernovae in the form- 
ing galaxies. While Ponman et al. (1999) investigate the 
entropy outside the core region, in §4 we investigate that 
in the central region of clusters where the information of 
cluster formation is expected to be kept although cooling 
is effective . 

The final issue in this paper is the metal abundance. 
The metal abundance of clusters is useful to know the 
chemical evolution of clusters and a possible correlation 
with dynamical quantities is a very interesting subject. 
ASCA has revealed the metal abundance distribution in 
clusters. Although the iron abundance outside the core 
region is roughly 1/3 solar and it is universal (Renzini 
1997; Fukazawa 1998), some clusters have abundance ex- 
cess (especially in iron) within the core region (Fukazawa 

1997) . However, not all clusters have the excess, and 
what makes the variation in the abundance excess is not 
well understood. Metals are ejected from galaxies in a 
cluster. Since the spatial distribution of galaxies in a 
cluster is more centrally concentrated than that of intr- 
acluster gas, and since a giant cD galaxy is situated at 
the center for most clusters, it is natural to expect that 
the metal abundance distribution has a central excess 
(Fujita, Kodama 1995; Ezawa et al. 1997). Thus, the 
non-existence of the abundance excess in some clusters 
may imply a mechanism to deform the original abun- 
dance distribution. However, the mechanism for this is 
still unknown. Fukazawa (1997) shows that the corre- 
lation between temperature of clusters and metal abun- 
dance at the centers is weak, while Allen and Fabian 
(1998) indicate that clusters with cooling flows tend to 
have the abundance excess. In §5, we investigate the 
relation between the metal abundance at the center of 
clusters and the virial density. 

Based on the above results, we discuss the evolution 
of the core region of clusters in §6. Section 7 is devoted 
to the conclusions. Throughout this paper, we assume 
H a = 50 kms^Mpc" 1 . 

2. Fundamental Plane Analysis with Einstein 
Data 

In paper I, we found that clusters form a fundamental 
plane and band in the (logp 0} \ogR 7 \ogT) space, using 
ROSAT data of 45 clusters with z ~ obtained by Mohr 
et al. (1999); the sample is nearly flux-limited. In this 
section, we show that the relations found in Paper I are 
also reproduced with Einstein data and prove that the 
relations, especially the one between po and R, are not 
significantly affected by observational errors such as un- 
certainties in the treatment of a central excess in emission 
and a correction of point spread function of telescopes. 
We use Einstein data of p and R for 46 X-ray clusters 
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with z ~ obtained by Jones and Forman (1984), al- 
though it is not a complete sample. Twelve clusters over- 
lap between the catalogue of Jones and Forman (1984) 
and that of Mohr et al. (1999). The densities and core 
radii are obtained by fitting surface brightness profiles by 
the conventional [3 model, 

Po 



[1 + (r/i?) 2 ]3/3/2 



(1) 



where r is the distance from the cluster center and f3 
is a fitting parameter. If a central excess in emission is 
seen, Jones and Forman (1984) excise the photon count 
data of the innermost region of the cluster, when they 
fit the surface brightness profiles to equation (|l|). Thus, 
the effect of local luminosity excess is excluded. On the 
other hand, Mohr et al. (1999) fit the excess with another 
[3 model instead of excising it, noting that in the way 
of Jones & Forman (1984), the best fit core radius is 
correlated with the size of the excised region and the 
result is rather subjective. But in this section, we simply 
follow the prescription of Jones & Forman (1984). 

The data of temperature are taken from the ASCA 
data catalogues if the clusters are observed by ASCA 
(17 clusters; Fukazawa et al. 1998; Markevitch 1998). 
Since not all clusters in the catalogue of Jones & Forman 
(1984) are observed by ASCA, the rest is taken from 
the catalogue of White, Jones, & Forman (1997), where 
they derived the temperatures using Einstein data by the 
deprojection method (Fabian et al. 1981). Although in 
paper I, we used the temperature data of previous ASCA, 
Ginga and Einstein observations gathered by Mohr et al. 
(1999), no significant difference in the results is found by 
this choice. 

As 

in Paper I, the data plotted in the (logpo, logi?, logT) 
space are fitted with a plane, 



Alogpo + B\ogR + C\ogT + D = 0. 



(2) 



The result of the least square fitting with equal weight 
for simplicity is A : B : C = 1 : 1.40 : -1.48. The scatter 
about the plane is 0.08 dex. In paper I, we obtained the 
result of A : B : C = 1 : 1.39 : -1.29 with a scatter 
about the plane of 0.06 dex. The scatters are nearly 
equal to typical observational uncertainties. Thus, the 
existence of the 'fundamental plane' is also confirmed for 
the Einstein data, and it is consistent with the ROSAT 
data. From now on, we sometimes represent Einstein and 
ROSAT results with indices E and R, respectively. 

In order to study the distribution of the observational 
data on the fundamental plane, we fit the data to another 
plane. 



a log po + b log R + c log T + d = . 
under the constraint, 



(3) 



This means that the plane (^) is perpendicular to the 
fundamental plane (Q). The result is oe : &e '■ ce = 1 : 
1.81 : 2.38. The scatter about the plane is 0.2 dex. In 
paper I, we found that or : 6r : cr = 1 : 1.18 : 2.04, with 
a scatter about the plane of 0.2 dex. We call this plane 
'the vertical plane'. Two unit vectors perpendicular to 
the fundamental and vertical planes are respectively de- 
fined by 



1 



ei 



e 2 



VA 2 + B 2 + C 2 



(A,B,C) 



1 



\Ja 2 + b 2 + c 2 



(a,b,c) 



(5) 



(6) 



For Einstein data, e±, and one of the unit vectors 
perpendicular to them are respectively given by 



eiE = (0.44,0.62,-0.65), 
e 2E = (0.32,0.57,0.76), 
e 3E = (0.84,-0.54,0.057). 
For comparison, the ROSAT data give 
em = (0.47,0.65,-0.60), 



e2R 



(0.39,0.46,0.80), 



e 3 R = (0.79,-0.61,-0.039) 



(7) 
(8) 
(9) 



(10) 

(11) 
(12) 



(Paper I). Since the observational uncertainty of each 
component of vectors is ~ 0.1, the sets of vectors (0)- 
(0) and ([Io|)-(|l2]) are consistent with each other. For 
Einstein data, the equations X E = p^ Ai R 0M2 T' - 65 , 
Y E = p Q 32 R°- 57 T°- 7e , and Z E = pg.MR-o.B4jiO.067 are 
three orthogonal quantities. Similarly, for ROSAT data, 
A R = pO- 47 i? - 65 T- - 60 , Y R = p° 39 R 0A6 T - 80 , and Z R = 
PQ 79 i? _0 - 61 T _0 039 are three orthogonal quantities. Fig- 
ure la shows the cross section of the fundamental plane 
viewed from the Ye axis, while figure lb shows the data 
on the (Ye j Z~e) plane. These figures correspond to figures 
1 and 2 in Paper I, respectively. As can be seen, a clear 
correlation ('fundamental band') exists on the plane. 

We compare the relations among physical quantities 
based on Einstein data with those based on ROSAT data. 
The vector e3E, which corresponds to the major axis of 
the fundamental band, means that 



Po oc R 



-1.6±0.3 



and 



rp „ n-0.1±0.2 „ 0.07±0.1 
1 OC ±1 OC /9q 



(13) 



(14) 



These relations should be compared to the results of 
ROSAT data obtained in Paper I: 



Aa + Bb + Cc = . 



(4) 



po oc R 



-1.3±0.2 



(15) 
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and 



T z?0.06±0.1 „ -0.05±0.1 

l oc it. oc p 



(16) 



Thus, for both data sets, the major axis of the funda- 
mental band is nearly parallel to the log po — log R plane, 
i.e., temperature varies very little along the fundamental 
band. 

In paper I, we showed that the observed Lx — T re- 
lation corresponds to a cross section of the fundamental 
plane. Moreover, we investigated the physical meaning 
of the observed Lx — T relation. In the rest of this sec- 
tion, we show that the fundamental plane derived from 
the Einstein data has the same physical meaning as that 
derived from the ROSAT data. We here assume that the 
scatter around the fundamental plane is due to observa- 
tional uncertainties alone, that is, A log (or A log X-r) 
is essentially zero. In this physical quantity cor- 

responding to each point on the fundamental plane can 
be represented by two parameters. For example, if we 
represent the gas mass fraction / with virial mass M v ; r 
and density p v ; r , the result turns out to be 



/ ex M^p~°- 



(17) 



for the ROSAT data (Paper I). On the contrary, the anal- 
ysis using the Einstein data results in 



f oc M°V : ° 



(18) 



In relations fll^ ) and (|l8|), we assume the usual scaling 
relations M v i r oc RT and p V ir oc R~ 2 T for dynamical 
equilibrium. Since R is the core radius, M v - lr and p v \ v 
respectively represent the virial mass and density of the 
core, strictly speaking. Both ( |l7j ) and (18) show that the 
observed Lx — T relation can be interpreted as a relative 
paucity of gas in small clusters. 

Since we have confirmed that the Einstein and ROSAT 
data give the same relationships, and since the spatial 
resolution of ROSAT is superior to that of Einstein, we 
use the ROSAT data in the following sections. 

3. The Central Excess Emission 

Up to now, we have concentrated on the outer com- 
ponents of surface brightness profile, but it is interest- 
ing to examine the nature of the inner components, or 
the central excess emissions when they exist. In this 
section, we extend the fundamental plane analysis to 
the inner components. In the catalogue of Mohr et al. 
(1999), the core radii i?; n and the central gas densities 
of the inner components po,in {R2 and po in their nota- 
tion) are given. Although the emissions from the inner 
components include those from several temperature com- 
ponents, in most cases they are dominated by the emis- 
sion from higher temperature components whose temper- 
atures nearly equal to the average temperatures of the 



clusters (Tamura et al. 1996; Xu et al. 1997, 1998; Ikebe 
et al. 1997). Thus, we here assume that emission from 
the inner components is totally from the higher temper- 
ature components. In fact, Mohr et al. (1999) assumed 
that the clusters are isothermal when they derive R m 
and po,in- We later discuss the uncertainties of this as- 
sumption. For the temperature of the inner and outer 
components, we use spatially averaged temperature T av 
presented by Mohr et al. (1999). In figure 2, we lap the 
data of the inner components over figures 1 and 2 in Pa- 
per I. The definitions of the axes are the same as those 
in Paper I. The data of inner components are situated on 
the fundamental plane. 

We consider more quantitatively the band distribution 
as we did in Paper I and ^2]. We fit the combined data 
including those of the inner components to equations (|^) 
and (||) . When the combined data of Paper I and those 
of inner components are fitted, the unit vectors perpen- 
dicular or parallel to the planes are 



en = (0.43,0.60,-0.68), 
e 2 i = (0.46,0.50,0.74), 
e 3 i = (0.78,-0.63,-0.058), 



(19) 
(20) 
(21) 



where index I represents the sample including the in- 
ner components. We may define X\, Y\, and Z\ from 
equations (p^|)-(|2l|) in the same way as Xe, Ye, and Ze- 
Taking account of the fact that observational uncertain- 
ties of each component of the vectors are typically ~ 0.1, 
equations (|l9|)-(^l|) are consistent with equations jlC|) 
(|l2|). Thus, the data of the inner components satisfy the 
relations of the fundamental plane and band we found 
in Paper I, that is, satisfy the same relations as those 
among physical properties of the outer components. Fig- 
ure 2 shows that the values of Ar and Yr of the in- 
ner components are not much different from those of the 
outer components. However, the values of Zr of the in- 
ner components are significantly larger than that of the 
outer components, because p v i T (oc R~ 2 T) and p of the 
inner components are larger. 

It is to be noted that there may be systematic errors in 
the data we used, because we have used the values of the 
gas density and core radius derived on the assumption 
that clusters are isothermal. Recently, the structure of 
the intracluster gas at the cluster centers is investigated 
in detail with ASCA for several clusters. When a cluster 
has two components in the surface brightness distribu- 
tion, the inner component often has two temperatures of 
>3 keV and ~ 1.5 keV (e.g. Centaurus cluster; Ikebe et 
al. 1999). The detailed observations show that the hot 
and cold components coexist in the core (e.g. Fukazawa 
1994; Ikebe et al. 1999). Moreover, spatially resolved 
temperature profiles obtained by ASCA indicate that the 
temperature of the hot components is nearly equal to 
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T av (Fukazawa 1994; Xu et al. 1998; Ikebe et al. 1999). 
Thus, the assumption of isothermality is correct when 
the emission from the inner component is dominated by 
the one from the hot component. In fact, for clusters 
such as A1060 (Tamura et al. 1996), AWM7 (Xu et al. 

1997) , Hydra-A (Ikebe et al. 1997), and A1795 (Xu et al. 

1998) , the emissions from the centers are dominated by 
those from the hot components. On the other hand, for 
some clusters, such as Centaurus, both emissions from 
the hot and cold components comparably contribute to 
the central emission excess, although the hot component 
dominates in volume. Ikebe et al. (1999) investigate the 
central region of the Centaurus cluster in detail; they 
determine po,in and i?j n of the Centaurus cluster by the 
analysis combining the ASCA and ROSAT data. They 
distinguish p ,in and i?; n for the hot components from 
those of the cold components in the central region. More- 
over, they take account of the inhomogcncity of metal 
abundance contrary to the analysis by Mohr et al. (1999) . 
They find that po.in ~ 0.02 cm~ 3 and R m = 46 kpc for 
the hot component, and that po,in ~ 0.06 cm" 3 , R m = 12 
kpc for the cold component. The data obtained by Mohr 
et al. (1999) are consistent with those for the cold com- 
ponent but not for the hot component. This indicates 
that there exist systematic errors of factor 4 in po,in and 
i?i n of the inner component in this cluster. (It should 
be noted that even if the data obtained by Ikebe et al. 
(1999) are used, the data of inner hot component of the 
Centaurus cluster are still on the fundamental plane as 
is shown in figure 2.) 

However, the Centaurus cluster seems to be the most 
extreme case. Among the clusters whose central regions 
are closely observed by ASCA (Centaurus, Virgo, A1060, 
AWM7, Hydra-A, A1795, and Fornax), the Centaurus 
cluster has the strongest abundance gradient. Moreover, 
only Centaurus and Virgo (Matsumoto et al. 1996) have 
the X-ray emission from the cold components comparable 
to the hot components. Furthermore, the central emis- 
sion of the Centaurus cluster is situated on the high 
end of the fundamental band (figure 2), while those of 
A1060 and others arc not. Since po of clusters at this 
end is large and the cooling time of the intracluster gas is 
short, the clusters are expected to have strong emissions 
from the cold components. Thus, we conclude that for 
most clusters we investigated, the hot components dom- 
inate the central excess emissions and the excesses are 
related to the fundamental band. Therefore, for most 
clusters, the systematic errors of the inner components 
should be smaller than those of the Centaurus cluster. 
Future observations will clarify this point. 

4. Entropy 

The entropy of intracluster gas is intimately related 
to the thermal history of the cluster. Virializaton of a 
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cluster at the formation and additional heating processes 
produce the entropy. Ponman et al. (1999) investigated 
the intracluster gas at a certain radius outside core re- 
gion where radiative cooling is not effective. They fix the 
radius at 0.1 times the present virial radius (roi), thus 
taking a nearly fixed value of p gas - They find that hot 
clusters satisfy the relation, 

S gas = log T + const , (22) 

which would be satisfied when gravitational collapse 
alone produces the entropy and when the profiles of gas 
density are similar, while cool (T < 4 keV) clusters does 
not. Moreover, they find that the entropy has a mini- 
mum value (~ 100 keV cm 2 ), which they call the 'en- 
tropy floor' for cool clusters. They suggested that the 
floor is a relic of the energetic galactic winds before the 
cluster formation. 

On the other hand, since the gas is affected by cool- 
ing in the core region of clusters, we expect that the 
entropy there provides information about not only heat- 
ing but also cooling. Following the study of Ponman et 
al. (1999), we first consider the relation between the en- 
tropy and temperature of gas at the center of clusters. 
The entropy of the gas at the cluster center is given by 

S gas - log(p- 2/3 T) + const , (23) 

where po is the central gas density derived by f3 fitting 
(Mohr et al. 1999). If a cluster has two components in the 
surface brightness distribution, we discriminate between 
po of the inner components and that of the outer com- 
ponents by index in and out, respectively. Since Mohr 
et al. (1999) do not present gas densities of the outer 
components po.out, we give them by 

T ° Ut jj ~ ) PO.in , (24) 
Jin^tout / 

where 7 ut and 7; n are the central surface brightnesses 
corresponding to the outer and the inner components, 
respectively, and R out is the core radius of the outer com- 
ponent. 

Figure 3a shows the relation between the entropy and 
temperature of gas at the center of clusters. Contrary to 
the result of Ponman et al. (1999), it has a large scatter. 
This is because not only T but also po varies; po takes a 
wide range of values in comparison with the gas density 
at roi- Most of the variation of p reflects that of p vlr 
(Paper I), that is, the formation redshift of the cluster. 
Thus, our method enables us to treat the relation be- 
tween the entropy and the cluster formation redshift in 
contrast of Ponman et al. (1999). In order to investigate 
the relation, we consider the relation between S g3S and 
Pvir- The density of p v ir is assumed to be proportional to 
i?~ 2 T, where R and T are the core radius and tempera- 
ture of the cluster, respectively, thus represented in the 
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units of keV Mpc~ 2 . If a cluster has two components in 
the surface brightness distribution, we use the core radii 
of the components i?j n and i? ou t, respectively, and we use 
the common temperature T av . The data are taken from 
Mohr et al. (1999). 

Figure 3b shows the relation between p vil and 5| as (= 
10 Sgas ). As can be seen, the variation of Sg as is related 
to that of pvir- Moreover, the correlation is tight. This 
may mean that in the central region of clusters, the in- 
ner component prescribes the structure and evolution of 
intracluster gas. We fit the points in Figure 3a and 3b to 
the following plane, 

Si log p vlr + S 2 log T + S 3 log 5* as + Si = . (25) 

The result is (Si, S 2 , S 3 ) = (-0.4,0.9,-1), that is, 



Sgas = log(p 



0.4±0.1^i0.9±0.3-i 
vir / 



const. 



(26) 



Equation ( pq ) is a linear transformation of the funda- 
mental plane including the inner components (eqs.]l9|]- 
plj]), because this is nothing but the relation among 
po, R, and T. Indeed, if we use the relation (p"7|), 
which is equivalent to the fundamental plane in the 
(logp , logi?, logT) space, we obtain 



'-'gas 0< - 



T 



T 



T 



rp3/5 

* — ■ ( 27 ) 



This is nearly equivalent to the equation (|26|). We can 
represent 5| as , p vlT , and T by X\, Yj, and Zj, which cor- 
respond to equations (p^)-(|2l]), 

o* „ y— 1.0y0.4 7— 0.6 
'-'gas A l x l ' 



gas 



I ! 



T oc X-^Yj^Z- - 1 



(28) 
(29) 
(30) 



Since the scatters of Xj, Y\, and Z\ are A log At = 0.07, 
Alogli = 0.2, and AlogZi = 0.6, respectively, the vari- 
ations of Sgag and p v i r are mainly due to that of Z\. 
Thus, figure 3b nearly corresponds to figure 2a or a side 
view of the fundamental plane, because the difference 
between Zn and Z\ is small. Note that the relatively 
large uncertainty in the power of T in equation ( |26| ) and 
the slight difference of T-dependence between equations 
( p6| ) and (p7|) come from the fact that the scatter around 
the plane AlogJfi, which corresponds to observational 
uncertainty, is emphasized by the linear transformation, 
because S*Jp-° A cx X^Y^Zf 01 . 

The line in figure 3b is the gas entropy achievable 



through gravitational collapse alone (S* 



oc p. 



-2/3 



for 



/ (relation [^7j , jl?]] ) , that is, the central gas fraction of 
poor clusters (generally having large p v i r ) is smaller than 
that of rich clusters. Although the gas entropy at the 
cluster centers is larger than achievable through gravita- 
tional collapse alone for large p V i r clusters, it does not 
converge to a floor value, contrary to the result of Pon- 
man et al. (1999). This means that the cooling is effec- 
tive in the central region of the clusters. 

5. Metal Abundance 

In this section, we study the relation between central 
metal abundance and p v i r - Since there is a giant cD 
galaxy at the center of most clusters, the metal abun- 
dance at the cluster center may mainly reflect the in- 
terstellar medium and the chemical evolution of the cD 
galaxy. We study 22 clusters overlapping in the cat- 
alogues of Mohr et al. (1999) and Fukazawa (1997). 
We investigate the iron abundance ^Fo,in and the sili- 
con abundance .Zsi,in for the central region (here sub- 
script 'in' donates the central value rather than the inner 
components). The data of metal abundance are taken 
from Fukazawa (1997). We treat the inner component if 
a cluster has two components in the surface brightness 
distribution, as well as the clusters for which a single 
component fit is acceptable. As we did in the previous 
sections, we assume that clusters are isothermal. How- 
ever, for clusters with relatively strong emission from the 
cold component at the cluster center, the data may suffer 
from systematic errors as discussed in §3. For example, 
the Centaurus cluster should have p vlr between the two 
data in figure 4a, because the core radius of the inner 
component of the cluster derived by Mohr et al. (1999) 
nearly equals to that of the cold component (see. §3). 

Figure 4a shows that i?Fe,in is an increasing function 
of p v i r , while ^si.in reveals no clear tendency. The cor- 
relation between Z-p e ^ n and p v j r indicates that our ar- 
guments about the fundamental relations among po, R, 
and T are physically meaningful. If the range of variation 
Pvir oc R~~ 2 T were due to an observational inaccuracy of 
R, the correlation between Z-p e ^ n and p v i r would not ex- 
ist. Although the data of silicon abundance have a large 
scatter (~ 0.5 solar, figure 4a), the silicon to iron abun- 
dance ratio ^si,in/-^Fc,in is a decreasing function of p v ir 
(figure 4b). Since iron and silicon are supposed to be 
supplied by different sources (Type la and II supernovae, 
respectively) , this correlation suggests that the contribu- 
tion of Type la SNe is larger for clusters with larger p v iD 
while that of type II SNe is more or less the same for all 
clusters. 



a given T). The normalization of the line is adjusted to 
the clusters with small p v i r , because the influences of ad- 6. Discussion 
ditional heating and cooling are expected to be small for 
them. The observed entropies of large p v i r clusters are lo- 
cated above the line. This is because of the non-constant 



In §3, we find that for clusters with a double surface 
brightness profile, the inner components satisfy the same 
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relations of the fundamental plane and band as those of 
the outer emission components. One possible interpre- 
tation is that the inner and outer emission components 
correspond to the separate dark matter components col- 
lapsed at different epochs under a hierarchical cluster- 
ing scenario, and that the temperature of the inner hot 
component represents the depth of the potential well of 
the inner dark matter component as suggested by recent 
ASCA observations. Assuming that the core of a clus- 
ter forms when a major merger (we simply say 'collapse', 
hereafter) occurs and that the envelope grows through 
minor mergers or accretion after that (e.g. Salvador-Sole 
et al. 1998), the interpretation implies that the collapse 
has occurred at least twice in the cluster with two com- 
ponents in the surface brightness distribution. In other 
words, the inner component of dark matter formed by 
the first collapse and had grown through accretion; when 
the second collapse occurs, the outer component of dark 
matter formed, and the cluster has grown further through 
accretion. In this scenario, the inner component of dark 
matter has survived the second collapse. In figure 5, 
we predict the formation epoch of each component of 
clusters using the spherical collapse model (Tomita 1969; 
Gunn, Gott 1972, see Paper II). As is seen, for a cluster 
with double distribution in the surface brightness, the 
inner component formed at z ~ 4 ahead of the forma- 
tion of the outer component regardless of cosmological 
parameters; the outer component forms at z <~ — 0.5. 
Since the gravitational mass is proportional to RT, fig- 
ure 5b shows that the gravitational mass of the inner 
component is about 1% of that of the outer component. 
This implies that the clusters had already grown up to 
1% of their present masses at z ~ 4. These may be con- 
sistent with the recent suggestion of Loken et al. (1999) 
and Miller et al. (1999) that clusters with central excess 
emissions reside in high density region of the Universe 
and started to collapse early. 

Figure 3b shows that the gas entropy at the center 
of large p v - lr clusters is larger than achievable through 
gravitational collapse alone (the line in figure 3b). On 
the other hand, figure 3b shows that, the gas entropy 
of large p v ir clusters is smaller than the 'entropy floor' 
(~ 100 kcV cm 2 ) claimed by Ponman et al. (1999). If 
the intracluster gas was heated by galactic processes, as 
a result of which the entropy floor was established be- 
fore the cluster formation as claimed by Ponman et al. 
(1999), our results suggest that the radiative cooling be- 
comes effective at the center of these clusters after they 
collapse, but that it is not so effective to cancel out the 
initial heating in spite of the short cooling time of the 
gas at the cluster center ( < 10 9 yr), which is well be- 
low the cluster age. Thus, heating sources should exist 
in the central region of the cluster even after the cluster 
formed. Considering the fact that the band distribution 
in figure 3b is tight and all kinds of components are lo- 



cated on the same band, the source must be stable and 
related to the evolution of the whole cluster. If the heat- 
ing process is temporary and local, the points in figure 3b 
would not form a band. 

From the above arguments, we may reject heating from 
the cosmic ray from AGN as the energy source, because 
the duration of the AGN activity is very short (~ 10 8 yr) 
and it cannot supply energy constantly. In fact, not all 
clusters have AGNs at the centers. Moreover the energy 
loss time of cosmic ray particles is well below the clus- 
ter age (e.g. Rephaeli, Silk 1995). Although supernovae 
seem to be more stable energy sources, the amount of 
energy is insufficient. For elliptical galaxies, which dom- 
inate in the central region of clusters, the supernova rate 
is at most 0.1 per 10 10 Lb© per 100 yr (Turatto et al. 
1994), where Lb© is the B-band solar luminosity. Since 
typical B-band luminosity of a cD galaxy is ~ 10 11 Lb©, 
the heating rate of supernovae is ~ 10 42 erg s _1 , assum- 
ing that each supernova ejects the energy of 10 51 erg. 
This is smaller than typical X-ray emission in the core 
region of a cluster (<~ 10 43 ~ 44 erg s _1 ). Multi-phase cool- 
ing flow (Fabian 1994) and heat conduction (Takahara, 
Takahara 1979) may stably transport energy form out- 
side region of a cluster, where high-entropy gas should 
be generated by shock heating through accretion (Me- 
zler, Evrard 1994; Navarro et al. 1995). In particular, 
the former model may explain the existence of the cold 
components in the central region of clusters. However, for 
the former model, the ultimate fate of the cooled gas is 
still unexplained (e.g. McNamara, Jaffe 1994; Voit, Don- 
ahue 1995). For the latter model, a mechanism of local 
reduction of heat conductivity is needed. This is because 
the classical conduction rate is too large to allow the exis- 
tence of the cold gas component (Sarazin 1986; Bregman, 
David 1988), while such a component certainly exists in 
the core regions of clusters (Ikcbe et al. 1999; Fukazawa 
1997). Therefore, there is no satisfactory model of the 
energy source or transmission so far. 

In §5, we investigate the relation between metal abun- 
dance (Zsi,in and ^Fe,in) and pvir. Since iron is mainly 
produced by Type la SNe and silicon is mainly produced 
by Type II SNe, this reveals the difference between the 
history of Type la and Type II SNe rates. Figure 4a 
shows that there is no correlation between p vlr and i?si,in- 
This may indicate that most of Type II SNe occurred be- 
fore the collapse of the central region of clusters. On the 
contrary, figure 4b shows that ^si,in/-^Fe,in is a decreas- 
ing function of p v ir- This decrease is chiefly attributed 
to the fact that the iron abundance excess is seen at the 
center of the clusters with large p v i r . This implies that 
the iron ejected from the galaxies at the cluster centers 
has accumulated in large p v ir and thus early collapsed 
clusters, although other possibilities cannot be ruled out. 
The clusters with p v - lr > 10 3 keV Mpc -2 especially have 
low silicon to iron abundance ratio (figure 4b). When 
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flo = 1 (pvii oc [1 + z] 3 ), this means that the central 
iron abundance excess is observed only in the clusters 
which collapsed at z > zp c ^2 — 3, assuming that the 
left end of the distribution (/9 vir <~ 20 keV Mpc~ 2 ) corre- 
sponds to the clusters which collapse at z ~ 0. When a 
cluster collapses, we expect that the existing iron excess 
at the center is smoothed out because of some violent 
mixing processes of intracluster gas. Thus, the obser- 
vational data may show that the iron had been ejected 
from Type la SNe into the intracluster gas at the clus- 
ter center for z > zp c , if the evolution of Type la SNe 
rate is universal. In other words, if a cluster collapsed at 
z = z co \\ > zpc, the iron excess had grown again during 
Zcoii > z > zpc, although it was once smoothed out at 
Zcoii- On the other hand, if clusters formed at z < zp c , 
little iron is supplied from Type la SNe after the forma- 
tion and thus no iron excess is observed. Assuming that 
the galaxies in the clusters formed at z >> z-p c , these 
imply that the iron ejection had continued for the ini- 
tial few Gyrs. The time-scale is roughly consistent with 
the theoretically predicted one (e.g. Mihara, Takahara 
1994). This is almost the same for ilo = 0.2. Finally, we 
comment on the relation between cooling flows and the 
central abundance excess suggested by Allen and Fabian 
(1998). Since old clusters have compact and dense cores, 
their central gas density also tends to be large. Thus, the 
clusters satisfy the criterion of cooling flow, that is, the 
cooling time of gas is smaller than the Hubble time. 

7. Conclusions 

We consider the scaling relations of the core region 
of clusters using the X-ray data. We first confirm that 
the fundamental relations among central gas density, core 
radius, and temperature we found using ROSAT data are 
reproduced by Einstein data, too. This means that the 
relations, especially between density and core radius, are 
not significantly affected by observational errors such as 
uncertainties in exclusion of the central emission excess 
and a correction of point spread function. This result 
strengthens our interpretation that clusters form a two 
parameter family in terms of mass and virial density. 

We find that the central excess in emission observed 
in many clusters also satisfies the fundamental relations. 
If the excess reflects the existence of the distinct dark 
matter component from the whole cluster as suggested 
by Ikebe et al. (1996), our result implies that the inner 
component collapsed ahead of the collapse of the whole 
cluster in the scenario of hierarchical clustering and that 
it represents the gravitational potential of the compo- 
nent. 

Gas entropies at the centers of old clusters show that 
both cooling and heating affect the thermal evolution of 
the intracluster gas. Moreover, the relation between the 
entropy and the virial density suggests that stable energy 
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sources or stable energy transmission mechanisms should 
exist in the central region of clusters. The ratio of sili- 
con to iron abundance at cluster centers is a decreasing 
function of the central virial density. This can be inter- 
preted that the iron ejection from Type la supernovae 
had occurred mainly before z ~ 2 — 3. 

The upcoming X-ray telescopes with high spatial and 
spectral resolution such as Candra, XMM, and ASTRO- 
E will be quite useful for detailing this kind of studies. 

This work was supported in part by the JSPS Research 
Fellowship for Young Scientists. 
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Fig. 1.. The observational data projected (a) on the Ze — Xe 
plane, (b) on the Ye — plane. 
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Fig. 2.. The observational data projected (a) on the Z-g_ — Xr 
plane, (b) on the Yr — Zjj plane. Filled squares, open squares, 
and crosses represent the inner components, the outer compo- 
nents, and the clusters without double distribution of surface 
brightness distribution, respectively. A diamond shows the hot 
inner component of the Centaurus cluster obtained by Ikebe et 
al. (1999). 
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Fig. 3.. The gas entropy at the cluster center as a function of 
(a) temperature of the cluster and (b) central virial density. 
Symbols are the same as those in figure 2. A typical error bar 
is shown at the upper right of the figures. A dashed line in 
figure 3b is proportional to p vi ^ 3 - 
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Fig. 4.. (a) The metal abundance at the cluster center as a 
function of central virial density of the cluster. In figure 4a, 
triangles and crosses represent Zp c i n and Zgi ; n , respectively. 
The point of iron abundance with an open circle represents the 
cluster with a cD galaxy (Fukazawa 1997). The uncertainties of 
•Zsi.in are typically 0.5 solar. A diamond shows the hot inner 
component of the Centaurus cluster obtained by Ikebe et al. 
(1999). (b) The abundance ratio ^si,in/^Fc,in is represented as 
a function of central virial density of the cluster. 



Fig. 5.. Theoretical predictions, (a) Radius-density relation (b) 
radius-temperature relation (c) density-temperature relation. 
Solid line: Q = 0.2 and M vir (t ) = 10 16 A%. Dotted line: 
Cl = 0.2 and M vir (t ) = 5 X 10 14 A^. Dashed line: U = 1.0 
and M vil (t ) = 10 16 ^g. Dash-dotted line: fi = 1.0 and 
M vil (t ) = 5xl0 14 M 9 . The open diamonds, triangles, and cir- 
cles correspond to the collapse redshifts of z co n = 0, z co n = 0.5, 
and z co n = 4, respectively. The observational data (po, R, and 
T) are overlaid being shifted moderately in the directions of po 
and R (/p v ir = 0.06po an< J r v - lr = 8R). Filled squares, open 
squares, and crosses represent the inner components, the outer 
components, and the clusters without double distribution of 
the surface brightness, respectively. The error bars of the last 
two are omitted but presented in Paper I. 



